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The possible occurrence of either a charge-density-wave 
or a Kohn anomaly is governed by the presence of Fermi- 
surface nesting and the subtle interaction of electrons and 
phonons. Recent experimental and theoretical investigations 
suggest such an effect for the hydrogen covered Mo and W 
(110) surfaces. Using density-functional theory we examine 
the electronic structure and the electron-phonon coupling of 
these systems. Besides good agreement with the experimen- 
tal phonon frequencies our study provides a characterization 
and quantitative analysis of an unusual scenario determining 
the electronic, vibrational, and structural properties of these 
surfaces. 

PACS numbers: 63.20.Kr, 68.35. Ja, 73.20.At, 73.20.Mf 

It is well established that the vibrational spectra of 
some transition-metal crystals exhibit more or less pro- 
nounced anomalous properties. For example, in the 
phonon band structure of bulk Mo a clear dip occurs at 
the 77-point Several theoretical studies |,|,|,|j7| 

support the idea that these anomalies are directly related 
to a nesting of the Fermi surface 

If the electron-phonon (e-ph) coupling and the Kohn 
anomaly become strong, an instability will arise, lead- 
ing to a Peierls-type, or more generally to a charge- 
density-wave (CDW) type of insulator. The new ground 
state then contains a frozen-in periodic lattice distortion. 
The distortion lowers the kinetic energy of electrons near 
the Fermi surface through an energy gap, which acts as 
the CDW order parameter. In systems of effectively re- 
duced dimensionality, such as layer and chain bulk com- 
pounds, and crystal surfaces, one expects nesting, and 
thus the occurrence of either strong Kohn anomalies, 
or of a proper CDW ground state, to be generally en- 
hanced Nonetheless, with the exception of W(IOO) 
and Mo(lOO) , no well-established cases of either giant 
Kohn anomalies or of surface CDW's have yet emerged 
on real crystal surfaces. 

The discovery of extremely sharp and deep indenta- 
tions in the surface phonon spectra of H/Mo(110) and 
H/W(110) at full hydrogen coverage has therefore 

come as a surprise, raising much interest, and many ques- 
tions. The anomalies are seen at an incommensurate k- 
point, Qci, along the [001] direction (TH) and, somewhat 
less visibly but also unmistakably, at the commensurate 



zone boundary point Qc2 = S of the surface Brillouin 
zone. Out of the ordinary Rayleigh mode, two simulta- 
neous anomalies appear to develop at these two k-points. 
One, uji, is extremely deep and sharp, and is only seen by 
helium atom scattering (HAS). The other, UJ2, is instead 
soft and mild, and is observed by both HAS and high res- 
olution electron energy loss spectroscopy (HREELS) p3| . 
The agreement between HAS and HREELS dispersion 
curves of the mild anomaly UJ2 supports the idea that 
this indentation is associated with the Rayleigh wave. 
The interpretation of the deep and sharp anomaly is less 
straightforward. A close link between the surface phonon 
anomalies and hydrogen vibrations jl3| seems to be ruled 
out, since the HAS spectra remain practically unchanged 
when deuterium is adsorbed instead of hydrogen ||ll|,|l2| . 

Recent density-functional-theory calculations of the 
Fermi surface of the hydrogen-covered Mo(llO) (1 x 1) 
surface identified a hydrogen adsorption induced quasi- 
one-dimensional Fermi-surface nesting in two different 
directions [|l4|,^. The calculated nesting vectors Q*i 
along TH and Q^J^ along TS are in excellent agreement 
with the critical wave vectors at which HAS detected the 
anomalous behavior. 

These experimental and theoretical results leave two 
possible scenarios still open for the ground state of these 
surfaces, namely: 

A. The surface e-ph coupling is very strong, the Ray- 
leigh-wave frequency at 2kF becomes imaginary, 
and the surface is unstable. A stabilizing distortion 
takes place and the true ground state of the sur- 
face is a quasi-one-dimensional CDW, incommen- 
surate along TH, and commensurate along TS. In 
this hypothesis, suggested in Ref. |l^, the surface 
phonon anomalies loi and 102 are therefore identified 
as phase and amplitude CDW modes, well studied 



in one-dimensional cases 



B. The e-ph coupling and the associated Kohn ano- 
maly UJ2 are only moderate. Therefore, the 
unreconstructed surface structure remains stable 
and the electronic instability is not removed. 
One should find only a slight downward shift of 
the Rayleigh-wave frequency which remains real. 
Within this picture, which was put forward in 
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Ref. [Q, the mode UJ2 is predominantly phonon- 
like and the other mode, uji, is predominantly due 
to an electron-hole excitation. 

Each of the two possibilities would be interesting and 
uncommon. In order to ascertain whether scenario A or 
scenario B is actually realized, a microscopic calculation 
of the relevant e-ph coupling, and of the ensuing shifts of 
surface phonon frequencies, is needed. A frozen-phonon 
calculation at Qci is out of question, due to incommen- 
surability. Luckily, point S is instead simply commensu- 
rate, and here the calculations are feasible. 

We performed density functional theory calculations 
using the local-density approximation for the exchange- 
correlation energy functional For the self- 
consistent solution of the Kohn-Sham equations we em- 
ployed the full-potential linearized augmented plane- 
wave method Q. Our code enables the direct calcula- 
tion of atomic forces and, with a damped Newton dy- 
namics scheme, the efficient determination of the relaxed 
atomic structure |Q. The substrates are modeled by 
five and seven layer slabs repeated periodically and sepa- 
rated by a vacuum region of thickness equivalent to four 
substrate layers. For the potential the (Z,m) representa- 
tion within each muffin tin (MT) sphere is taken up to 
^max — 3 while the kinetic-energy cutoff for the interstitial 
region is set to 64 Ry. We choose the plane- wave cutoff for 
the wavefunctions to be 12 Ry and employ a {l,m) rep- 
resentation in the MTs with /max — 8 for them. The MT 
radii for the W and Mo atoms are chosen to be 1.27 A. For 
the hydrogen this radius is set to 0.48 A. Fermi smearing 
with a broadening of kT^i = 0.07 eV is used in order to 
stabilize self-consistency and k summation [p3[. In the 
case of the W surfaces the valence and semi-core (core) 
electrons are treated scalar (fully) relativistically while 
for Mo all electrons are treated non-relativistically. The 
in-plane lattice constants 3.14A (W) and 3.13A (Mo) 
calculated without including zero-point vibrations are in 
good agreement with the respective measured bulk lattice 
parameters (3.163 A and 3.148 A for W || and Mo [||, 
respectively) at room temperature. 

The necessary first step of our study is to determine 
the atomic and electronic structures of the clean and 
hydrogen-covered (110) surfaces. In Table | we summa- 
rize the calculated relaxation parameters for these sys- 
tems. It turns out that the results for W(llO) are very 
similar to those obtained for Mo(llO). The investiga- 
tions shine light on a long-standing problem concerning 
the structure of W(llO) covered with more than half a 
monolayer of hydrogen. For this system Estrup's group 
observed a symmetry change in the low energy elec- 
tron diffraction pattern and proposed that this might be 
caused by a displacement of the top layer W atoms along 
the [110] direction Similar studies for H/Mo(110) 

did not provide any evidence for such a top-layer-shift 
reconstruction p^. For both adsorbate systems we per- 



formed several structure optimizations starting from dif- 
ferent trial configurations in order to check whether such 
a reconstruction is energetically favorable. It turns out 
that for both substrates the hydrogen atom relaxes into 
a quasi-threefold position (indicated as H in Fig. 1 in 
Ref. ||l^,|l^), and that the adsorption reduces the inward 
relaxation of the clean surface considerably. Further- 
more, for H/Mo(110) |l|,|l5| as well as H/W(110) |2|] 
we find no evidence for a pronounced top-layer-shift re- 
construction. However, yi is non-zero, but very small, for 
both adsorbate systems, which indicates an anisotropic 
vibration of the top layer parallel to the substrate with 
a favorable direction along [TlO]. 

The similarities between Mo(llO) and W(llO) con- 
tinue when we compare the electronic structure of the 
relaxed surface systems. For both systems the H adsorp- 
tion alters the surface potential and induces a shift of the 
(dsz^-r^jdxz) band to higher binding energies |l^ , [l5| , p8| . 
In this way the Fermi line associated with this band is 
moved away from TN into a band gap of the surface 
projected band structure, and the respective states be- 
come true surface states. For both adsorbate systems 
the shifted {d^^2_j,2^ dxz) band is characterized by a high 
density of states at the Fermi level, and the new Fermi 
contour gives rise to a nesting . The magnitudes 

of the calculated and measured critical wavevectors along 
TS and TH are listed in Table |ll[ The agreement between 
theory and experiment for both systems is excellent. 

The experimental phonon spectra of the (110) surfaces 
of Mo and W show two distinct reactions to the ad- 
sorption of hydrogen. Along TH and TS a monolayer 
of hydrogen induces a softening of the Rayleigh and (to 
a smaller amount) of the longitudinal |2^ wave while a 
stiffening of these modes is observed along TN. Our goal 
is to investigate both effects theoretically. 

A frozen phonon calculation at S is particularly con- 
venient, since at the zone boundary the Rayleigh- wave 
polarization is purely vertical, and the second layer is 
immobile by symmetry. In fact, at S the Rayleigh wave 
is not embedded in the projected bulk phonon bands and 
thus it is strongly localized at the surface. Moreover, we 
find that the hydrogen vibrations are of no significant 
importance for the calculated phonon frequencies. 

In order to make our frozen phonon studies compu- 
tationally feasible we use an enlarged surface unit cell 
together with a five layer slab and reduce the plane- wave 
cutoff for the wave functions to 10 Ry (the convergence 
of total energy differences was tested). The geometries 
are defined by the five layer relaxation parameters pre- 
sented in Table ^. One should mention that some relax- 
ation parameters, i.e., Ad^, change considerably when 
we perform a calculation with a slab of five instead of 
seven metal layers. We note, however, that the values 
of the critical wavevectors and the position of the hydro- 
gen with respect to the substrate surface are practically 
insensitive to the thickness of the slab. This indicates 
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that the physical properties we are interested in, e.g. the 
nesting features, are well localized surface phenomena 
and that the results of our five layer slab studies can be 
trusted. 

The phonon frequency was evaluated using the har- 
monic term of a fourth-order polynomial fit to the total 
energy changes. The total energies were calculated for 
five displacements between 0% and 3.6% of the (110) 
surface interlayer spacing. One might expect that the 
calculated frequencies, especially the one for the S'-point 
phonon, are sensitive to the Brillouin-zone sampling. For 
instance, for the study of the (2x1) reconstructions of 
the diamond (111) surface Vandcrbilt and Louie |3^] em- 
ployed a k-point mesh which becomes logarithmically 
denser close to the zone boundary. We performed test 
calculations using 16, 56, and 120 k-points in the sur- 
face Brillouin zone in order to test the convergence of 
our results, and found that a uniform k-point mesh of 56 
points is sufficient. The calculated frequencies are col- 
lected in Table [11 . At the N point our results reproduce 
the experimentally observed stiffening of the Rayleigh- 
wave frequency as hydrogen is adsorbed. This agreement 
encourages to go on to the study of the S'-point Rayleigh 
wave. Here we find that the stiffening effects are over- 
compensated because of a strong interaction between the 
phonons and the electron states at the Fermi level. In 
our calculations this leads to a moderate lowering of the 
phonon energy by 5.5 meV for Mo and 6.1 meV for W, 
in good agreement with the experimental results. This 
strongly supports the conclusion that scenario B is the 
correct one. 

We may extract the surface e-ph coupling at the S 
point by calculating the splitting of the two Kohn-Sham 
eigenvalues folded back at the point midway between F 
and S" as a function of the nuclear displacement. We ob- 
tain the gain in electronic energy which arises from the 
lattice distortion of the frozen phonon, and this quantity 
can be related to the electron-phonon interaction at the 
surface. The calculated splitting 2Ae is practically linear 
in the nuclear displacement and the coefficient of the lin- 
ear relation is 2.81 ± 0.10 eV/A for Mo and 4.41 ± 0.06 
eV/A for W. The larger value for W reflects a more ef- 
ficient coupling between the electrons and phonons and 
suggests that the anomaly in the phonon band of W(llO) 
should be more pronounced. This is indeed confirmed by 
the HAS experiments, which find that the deeper dip 
is wider on H/W(110) than on H/Mo(110). Further- 
more, Hulpke and Liidecke ||ll|,|l^ reported the occur- 
rence of small structures in the HAS diffraction pattern 
for H/W(110) at the critical wavevector. These features, 
which were not observed for the Mo substrate, might in- 
dicate the start up of the formation of a CDW at W(llO). 

Our calculations of the electron-phonon interaction at 
the (110) surfaces of both Mo and W and their changes 
due to hydrogen adsorption pinpoint the phonon char- 
acter of the small anomaly uj2 and identify the subtle 



interplay between the electronic structure and the vi- 
brational spectra of the transition metal surfaces. The 
results clearly support the interpretation that the small 
dip observed by HAS and HREELS is due to a Kohn 
anomaly, while the huge dip arises from an electron-hole 
excitation. Thus, scenario B is operative. There are 
in fact two independent experimental hints which sup- 
port our conclusion. If picture A were correct, and uji 
a phason, that phason should exist at an incommensu- 
rate point like Qci, but should be absent at the zone- 
boundary commensurate point Qc2 = S However, 
data near S clearly indicate a lower branch at that point 
too. Only if scenario B holds, is a second branch ex- 
pected at S due to an electron-hole excitation and a good 
Fermi-surface nesting. Furthermore, recent HREELS ex- 
periments [ p9| reveal a positive temperature dependence 
of the small indentation {duj2/dT > 0) consistent with 
that of a Kohn anomaly, i.e., a weakening of the anomaly 
with increasing temperature. 

In conclusion, we have discussed the atomic and elec- 
tronic structures of the clean and H-covered Mo(llO) and 
W(llO) surfaces and studied the changes in the vibra- 
tional properties of these surface systems upon hydro- 
gen adsorption. Two decisive experimental findings could 
be reproduced: At the symmetry point N the Rayleigh 
phonon mode is stiffened by the hydrogen adsorption, 
whereas for the 5-point Rayleigh wave we obtain a weak- 
ening of the phonon energy due to Fermi-surface nesting. 
The calculations support the interpretation jl^ that sce- 
nario B is in effect (the Kohn anomaly picture), as op- 
posed to the CDW picture |jl^ (scenario A). Thus, the 
modest anomaly seen in HAS and HREELS is a predomi- 
nantly phonon-like excitation, and the giant indentation, 
seen only in HAS, is predominantly due to an excitation 
of electron-hole pairs. The various elements identified 
in our analysis could be combined in a future study of 
the behavior of the H adatoms in order to analyze their 
liquid-like state features observed in the HREELS mea- 
surements [|l3| . 
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TABLE I. Calculated relaxation parameters for the clean 
and H-covered (110) surfaces of Mo and W. The height of 
the hydrogen above the surface and its [TlO] offset from the 
[001] bridge position are denoted by dn and yu, respectively. 
The shift of the surface layer with respect to the substrate is 
yi. The parameters Adij describe the percentage change of 
the interlayer distance between the i-th and the j-th substrate 
layers with respect to the bulk interlayer spacing do . For each 
system the results for a five layer (first line) and for a seven 
layer (second line) slab are presented. 



system 


yu dn 

(A) (A) 


(A) 


Adi2 
(%do) 


Ad23 Ad34 

(%do) (%do) 


Mo(llO) 






-5.9 
-4.5 


-0.8 

+0.5 0.0 


H/Mo(110) 


0.63 1.08 
0.60 1.07 


0.05 
0.03 


-2.7 
-2.1 


-0.4 

+0.1 -0.1 


W(llO) 






-4.1 
-3.3 


-0.2 

-0.2 -0.6 


H/W(110) 


0.68 1.21 
0.67 1.20 


0.01 
0.05 


-1.3 
-1.3 


0.0 

+0.3 +0.4 


TABLE H. Theoretical Fermi surface nesting vectors com- 
pared to critical wavevectors obtained by HAS and HREELS 
experiments [pT|,p^jl^] . 


direction 


system 




IQc 










theory 


experiment 


FH 


H/Mo(110) 
H/W(110) 




0.86 
0.96 


0.90 
0.95 


rs 


H/Mo(110) 
H/W(110) 




1.23 
1.22 


1.22 
1.22 



TABLE HL Comparison of calculated frozen phonon ener- 
gies and experimental results values obtained by HAS [ pT|Jl^ 
and HREELS [^ . The numerical accuracy of the theoretical 
phonon energies is about ±lmcV. 



phonon 


system 


Eph 


(meV) 






theory 


experiment 


iV 


W(llO) 


17.0 


14.5 




H/W(110) 


19.8 


~19 


S 


Mo(llO) 


23.7 


~21 




H/Mo(110) 


18.2 


<16 




W(llO) 


19.1 


16.1 




H/W(110) 


13.0 


11 
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